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Abstract 

We perform the investigations of the resonant tunnehng via impurities embedded in the AlAs 
barrier of a single GaAs/AlGaAs heterostructure. In the liV) characteristics measured at 30mK, 
the contribution of individual donors is resolved and the fingerprints of phonon assistance in the 
tunneling process are seen. The latter is confirmed by detailed analysis of the tunneling rates and 
the modeling of the resonant tunneling contribution to the current. Moreover, fluctuations of the 
local structure of the DOS (LDOS) and Fermi edge singularities are observed. 
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I. INTRODUCTION 



A considerable amount of resonant tunneling studies was performed on single barrier 
GaAs/AlAs/GaAs heterostructures. AlAs is an indirect gap material with the minimum of 
the conduction band at the X point of the Brillouin zone, while in GaAs the minimum is at 
the r point. As one can see in Fig. 1, the bands are aligned in a way that AlAs layer forms 
a barrier for the electrons in the F valley (solid line), and a quantum well (QW) for the X 
valley electrons (dashed line). There were several analyses of tunneling through X- valley 
states including first observation of negative differential resistance in such structures ll, 
c„„siderat,o,. on T-X t.ansfe. .nechan... Q. d,scus.o„ of conse™t,o„ B a. 

well as through X-minimum-related donor states (including investigations of the splitting of 
the ground state 0,0], pressure coefficients 6] and binding energies j?! of donors). However, 
in those experiments macroscopic samples were used where many donors were involved. 
Recently, interest was more focussed on the observation of the tunneling through individual 
objects such as quantum dots or impurities. In this case it is possible to perform the 
spectroscopy of the two dimensional electron gas (2DEG) in the so-called emitter as well 
as of these individual objects which are involved in the tunneling, e.g. quantum dots P|. 

In this paper we report results of resonant tunneling experiments performed on 
GaAs/AlAs/GaAs single barrier junctions in which single impurities, incorporated in the 
AlAs layer, are involved. Some additional features (like Fermi edge singularity and fluctua- 
tions of the LDOS) are also discussed. 



II. SAMPLE FABRICATION 



Samples were grown by molecular beam epitaxy on a (lOO)-oriented Si-doped n-type 
GaAs wafer (n^ = 1 x 10^^ cm~^). The active part of the heterostructure consists of a 
10.2nm thick AlAs barrier incorporating a Si 5-doping in the middle with a concentration 
of = 1 X 10^° cm~^. In some reference samples the Si (5-doping is omitted. The barrier is 
separated from the heavily doped contacts {rid = 4 x lO^^cm"'^) by 200nm GaAs spacers in 
order to obtain high quality 2DEG. Mesa structures of different lateral sizes ranging from 
500 /xm down to 100 nm were fabricated. This allows us to have different numbers of donors 
within the junction, starting from more than 10^ for the larger mesas to less than 10 for the 
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smallest ones. 



III. EXPERIMENTAL RESULTS FOR LARGE MESAS 

Figure 2 shows current versus voltage characteristics measured at 4.2K on 500/im mesas, 
one (5-doped (10^°cm~^, solid line), the other one undoped (dashed line). In the liV) curve of 
the doped mesa we observe an important increase of the tunneling current at approximately 
1 Volt. This is related to the resonant tunneling through donor states j^, and since in 
this case the number of donors is huge (5 x 10^), the observed increase is very broad. The 
other current bump at about 2V, present for both samples, is due to the tunneling through 
X-minimum quantum well states 1^,1^]. 

IV. EXPERIMENTAL RESULTS FOR SMALL MESAS 

If the number of Si impurities inside the barrier decreases Q| (small mesas), then the 
contribution of one individual Si impurity can be resolved In Fig. 3 we present the I{V) 
characteristic of a 400nm mesa, measured at very low temperature (30mK). The diameter is 
so small that the expected average number of donors within the mesa is of the order of 10. 

We observe several sharp current increases (labeled for clarity), which are followed either 
by the decrease of tunneling current {e.g. 1,3) or by further current increases (2,4). The 
observed features can be qualitatively explained in the following way: each of the increases 
indicates the beginning of the resonant tunneling process via an impurity level, whose energy 
(e^) becomes aligned with the Fermi energy (ep) in the emitter. Then, as the voltage 
increases further, the impurity level goes down with respect to ep but remains aligned with 
occupied states in the emitter, so that the tunneling current is maintained. Once it is below 
the bottom of the subband in the emitter, the current decreases. 

The fact that we observe many steps spread in a relatively wide bias range indicates that 
some interdiffusion of Si atoms takes place during the growth and the Si 5-layer is broadened. 
The steps which appear at low biases correspond to donors placed close to the collector side 
of AlAs layer. In electric field, the position at which the binding energy of the impurity has 
its maximum is shifted from the centre of the AlAs layer [12, Il3| to the collector side of the 

n y ^ i 

structure |14| . 
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The above considerations lead to the conclusion that one should expect a rectangular 
shape in I{V) characteristics while measuring the tunneling from 2DEG through a localized 
state at very low T. However, a more complex structure is systematically observed. Let 
us discuss this point by taking as a reference step 3 (Fig. 0)). At approximately 825 mV 
a steep current increase indicates the beginning of the resonant tunneling, i.e. = tp. 
At the current onset, a sharp peak with an overshoot is observed. For higher values of 
bias (830-970 mV), the current rises systematically but several reproducible peaks with very 
small amplitudes are superimposed. Then at approximately 970 mV the tunneling process 
is switched off. At this bias = ei where ei is the bottom of the subband in the 2DEG (see 
Fig. 1). 

In order to confirm that the maximum of the current observed at the end of the step 
(at 970 mV) is related to the switching off of the tunneling and thus corresponds to the 
case where the impurity level is aligned with the bottom of the 2DEG subband in the 
emitter, we further performed magnetotransport measurements. I{V) characteristics have 
been measured with a magnetic field applied along the direction of the current. In Fig. El 
I{V) curves between OT and 2.4T are plotted and the shift of the peak related to the second 
Landau level is clearly seen as a function of the magnetic field (the first LL is hidden in the 
maximum at the end of the step, but it is revealed in high magnetic fields). When B — >■ 
OT they converge to ~ 970mV, which is a clear proof that this voltage corresponds to 
the bottom of the subband. It means that the total width of the emitter corresponds to 
140±10 mV. Using the correspondence between the Landau level separation measured in 
/(V") characteristics and the one calculated in energy scale huJc (where uJc is the cyclotron 
frequency), we were able to recalculate the width of the emitter into energy, obtaining ep-ei 
^ 4.8meV. 

Referring to Fig. we would like to note that within each I{V) curve a rich structure of 
very sharp peaks appears in magnetic field. The average separation between these peaks is 
of the order of O.lmeV and they are very sensitive to the magnetic field. Similar observations 
have already been reported in the literature ^. The origin of those structures is not fully 
understood yet and some further investigations must be performed. 

When the tunneling process stops, for voltages higher than 970 mV, a smooth decrease 
of the tunneling current is observed. It should be noted that the current never falls to its 
value observed before the beginning of the tunneling, but remains at a higher level. We 
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propose that this feature (current tail), together with the systematic rise of the current as a 
function of bias during the tunnehng, is due to the couphng to the environment (contribution 
of non- resonant tunnehng is neghgible in this scale), which is realized by interactions with 



phonons 1^. This issue will be discussed in the next paragraphs, where we present our 
model and results of the calculations. In the last part of the paper we will focus on the fine 
structures observed within the steps and discuss the current maximum at Fermi energy. 

V. MODEL 

Tunneling transfer times have been calculated taking into account only the interactions 
with the GaAs longitudinal acoustic phonons. Indeed, optical phonons have much higher 
energy (> 30 meV) and cannot be responsible for the observed features that appear at typical 
energies of only a few meV. To compute the tunneling time Tk, from a 2DEG eigenstate |kj > 
to the impurity we use the Fermi golden rule 

^ 271" ^ 

— = — |< v^l^oi-phlkj >f5{ef + hxj{q) - e^) (1) 

q 

where kj is the wavevector of the electron in the 2DEG plane, ej(j) is the energy of the initial 
(final) electron state and the kets \ip > and |kj > refer to the eigenstates of the impurity 
and of the 2DEG, respectively. In this equation the summation is performed over all the 
final states with one phonon of energy ^(q). Because our experiments are performed at 
very low temperature only phonon emission processes are relevant. 

The computation of the matrix elements requires the knowledge of the wavef unctions. 
We adopt the envelope function formalism in the single parabolic band approximation. 

The envelope eigenfunctions of the 2DEG are approximated by the so-called modified 



Fang-Howard wavefunction which is given by 
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<z,r\k>=ij{z)y'^-^ , (2) 

where r and z are the projections of the electron position vector in the 2DEG plane and in 
the direction of the current respectively, L is the lateral length and ip{z) is defined by 

ij{z) = Nzqc^'^"''^ if ^ < (AlAs barrier), (3) 
tl){z) = N{z + Zo)e"^"/^ a z>0 (GaAs). 



Here N, zq and b are defined by 



487re^me 
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^l + bzo + \bHl{l 



_b_ 



Zq = 



b + Kb^ 



h 



(4) 
(5) 

(6) 
(7) 



where mg and mi, are the electron masses in the well and in the barrier respectively: rUe = 
0.067mo, nib = 0.124mo, k is the dielectric constant in the well: k = 12.9 x kq, Ug is the 
electron density of states: Us — 1.32 x lO^^m"^ here, ridepi is the concentration of residual 
acceptors with a typical value of 6 x lO^^m"^, mo is the free electron mass and kq is the 
dielectric constant in the vacuum, is the barrier height and is the energy of the electron. 

For the sake of simplicity the wavefunction of the impurity inside the AlAs quantum well 
is approximated by a spherical gaussian function: 



< z,r\cp >^ cp{x)cp{y)cp{z - Zi) , 



where 



(fi{x) 



e 2^ 



(8) 



(9) 



Zi is the position of the impurity along the growth direction and a corresponds to the spatial 
extent of the donor wavefunction and is comparable to the Bohr radius of a Si impurity. In 
the AlAs barrier this value is estimated to be close to 2.6nm. The adopted form of the donor 
wavefunction is a very crude approximation, which cannot give good absolute values for the 
tunneling rates. However, we believe that it is good enough for the qualitative description 
of the physical phenomena involved. 

The electron-phonon interaction is given by: 



(10) 



where qy is the phonon wavevector in plane of the 2DEG, 6^ is the creation operator for a 
phonon of wavevector q (in all the space directions), and a{q) are given either by: 

1 



a{q) 



ee 



pz 



(11) 



for the piezoelectric scattering or by: 



for the deformation potential scattering, where D is the deformation potential constant, p is 
the GaAs density, Cg is the sound velocity in GaAs, 6*^ is the piezoelectric constant (13eV, 
5300kg/m^, 3700m/s and 0.14eV respectively), VL is the volume of the crystal. We have 
assumed that the phonons obey the Debye's law. 

The total phonon-assisted tunneling rate r is obtained by summing over the initial states 
|kj >: = J2k '''k^ usual calculation [3] gives: 



r h{27T 



dkiki I dO I dq\\q\\\f{qz)\^\a{q)\' 



ko 



hcs qz 



F\ki,q\\,e) , (13) 



where = ^y2meVa.aix{ed — ei, 0)/^, kp is the 2DEG wavevector at the Fermi energy ep, qo 
is the Debye wavevector, 6 is the angle between kj and qy, F is given by: 

F = 2y/naexp ~ ^^il^^ + kiq\\a'^ cos 9 

and f{qz) is a form factor defined by 



(14) 



f{q.) = / ipiz)^ljiz)e-"i^'dz (15) 

J — oo 

where integration can be performed over the whole space since cp and ip are localized. Finally 
qz is given by the energy conservation: 



and q'i + qfi = q^- 



VI. RESULTS OF THE CALCULATIONS 

In this section we present the numerical results. Figure El shows the computed phonon- 
assisted tunneling rates for both deformation potential and piezoelectric interactions. The 
rates are plotted as a function of the donor energy e^. The energy reference is the energy ei 
of the bottom of the 2DEG subband. Both interactions have the same qualitative behavior, 
but the deformation potential interaction gives a rate one order of magnitude higher than 
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that arising from the piezoelectric interaction. At the beginning, the rates increase with 
to reach a maximum close to ed=€i. This increase is mainly due to the fact that the number 
of electrons available for phonon-assisted tunneling increases when decreases. Finally, 
when ed diminishes further and goes below the 2DEG subband, the rates diminish too: the 
number of initial states available for tunneling remains constant but the tunneling rates with 
phonons of high energy are strongly reduced by the form factor f{qz) appearing in Eq. ^1 
This decreasing tail can be ascribed to the fact that only acoustic phonons with small q 
contribute strongly to the relaxation time in a 2DEG. 

In order to describe the observed shape of the tunneling current, one should take into 
account not only the phonon-assisted tunneling but also the elastic one. In the sequential 
model of the tunneling, one should consider two processes: transfer from the emitter to the 
impurity and from the impurity to the collector. In our system, the donors observed at low 
biases are located on the collector side of the barrier (even if the change of donor binding 
ene.gy at the edge of QW . ta.e„ into aeeount QQ). The total cu„e„t can be expressed 
by the tunneling rates between the emitter and the impurity only: 




where r^f and stand for the deformation potential and the piezoelectric potential contri- 
butions, respectively, r stands for the elastic tunneling contribution and is inversely propor- 
tional to the square of the overlap between the donor wave function and the electron wave 
function in the emitter. 

The amplitude of elastic contribution is adjusted to fit the experimental value of the cur- 
rent measured just above the threshold, where the contribution of phonon-assisted currents 
is still negligible and singularity ai ep already vanishes. If we take a donor placed close to 
the collector side of the barrier, e.g. Zi =7-9nm, we get a reasonable agreement with the 
experiment for a= 2.2nm-3.2 nm, respectively. The values of a are of the same order of 
magnitude as the Si Bohr radius. The result of the calculation for Zj=7nm and cr= 2.2nm is 
presented in Fig. [7| with a solid line. The calculated current rises up to = ei, in agreement 
with the experiment and the current tails at higher bias are reproduced. The calculated 
sharp decrease of the current at = ei is related to the step-like DOS in the emitter which 
is refiected in the elastic contribution. 

In order to improve the agreement with the experiment one has to find the reason why 
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the elastic contribution does not change so abruptly at = ei. This can happen due 
to e.g. the disorder-related broadening of the DOS. Indeed, in the magnetic field results 
(see Fig. Ej) - Landau level N=0 is always considerably broadened, almost independently of 
magnetic field. Possibly there is a band tail at the bottom of the conduction band induced 
by disorder. We can take it into account by including in our model the convolution of the 
unperturbed DOS with a gaussian. The result is shown in Fig. [7| with a dashed line. The 
FWHM of the gaussian was taken to be O.SmeV. One can also expect that more exact many- 
body calculations of the current would lead to such a broadening [l^l- Indeed, the Fermi 
Golden rule is an approximation that does not take into account tunneling events with more 
than one phonon. 



VII. THE FINE STRUCTURES 



Finally, we want to comment on the detailed structures observed on the steps, which 
have not been discussed until now. The first striking feature is the sharp maximum at the 
current onset, i.e. at the Fermi energy, see Fig. |3J This maximum is due to the Coulomb 
interaction between electrons in the 2DEG and the impurity state. This effect, known in 
X-ray spectroscopy, was first predicted for the case of tunneling through a localized level 
by Matveev and Larkin First experimental confirmations were presented by Geim 

et al. who performed careful temperature analysis in order to demonstrate a singular 
behavior of the current at the Fermi energy. In our case the system changes its properties 
due to charging effects during voltage sweep. Instead of overcoming these effects, we take 
advantage of them. In the I{y) curves a hysteresis of the tunneling current is observed: by 
sweeping up and down the bias, the same structures are observed in different bias ranges and 
they have different widths in the voltage scale. However, by applying a linear transformation 
to the bias and a vertical shift to the current (the need for the latter is related to charging 
of some capacitance during the voltage sweep), it is possible to superimpose the two liV) 
curves for both directions of bias sweep. This is shown in Fig. |H1 and Fig. El where the 
different axes are labeled for both bias directions. One clearly sees that the main difference 
is a shift of the current onset at the Fermi level in Fig. IHl It means that the number of 
electrons in the emitter was slightly modified. But no matter what the position of the onset 
was, the observed maximum remains stuck to the Fermi level. 
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/iner Q and Q 



FoUowine; [181 and [19|, the Fermi edge singularity has a power-law shape of the form 

/ oc (v/(eF-e)2 + r2) ~^ X + arctan Oiep - e) , (18) 

where e is the energy of the tunneling electron, h/Tc is the time an electron can stay in the 
impurity state before escaping into the collector side and 

/3 = |- {k^dy' , (19) 

where d is the average distance between the 2DEG and the impurity state and kp is the 
wavevector at the Fermi energy. In Fig. ^1 we present an enlargement of step 3 from Fig. 121 
for two polarizations. Despite the shift at lower biases, the two curves start to be very 
similar at biases higher than 0.85 V (in the lower scale of Fig. irUj) . It means that at biases 
higher than 0.85V the many-body contribution becomes negligible and the measured current 
is due to noninteracting electrons only. The amplitude of this current does not depend on 
the bias sweep direction. This limits the range of acceptable values of fitting parameters. 
The voltage corresponding to ep was taken from the experimental results and the value of 
kp = 9.3 X 10^m~^ was calculated from the measured ep-ei. Three parameters were fitted, 
as in Ref. the amplitude of the many-body current in Eq. , F^ and (3. The best fit 
was obtained for Fc ~ 0.01 meV and /5 ~ 1. The average distance d evaluated from (19) 
is close to 2.5 nm, which is only about 4 times smaller than expected for impurity placed 
at the collector side of the barrier and gives a reasonable order of magnitude. The transfer 
time T is of the order of lOOps. 

Another interesting feature that is worth pointing out is the existence of the hardly 
resolved peaks within the plateau, in Fig. El between 0.83V up to 0.94V (in the lower scale) . 



20 



We attribute them to local fluctuations of the DOS in the vicinity of an impurity 
Naturally, we can expect that if the fluctuations are different for some steps, then these steps 
are most probably related to different impurities and do not correspond to ground or excited 
states belonging to the same impurity. If we compare step 1 (Fig. ^ and step 3 (Fig. Ej), 
we see that the detailed characteristics are different, even if the signal is much more noisy 
in step 1, as the values of current are very small. Thus we argue that we do really observe 
fluctuations of LDOS and that step 1 and 3 correspond to two different impurities. 
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VIII. CONCLUSION 



We have studied resonant tunneling from a 2DEG through the single impurity state. We 
see contributions from individual impurities - each of them appears to be different from 
the expected rectangular shape. In particular, there is a current enhancement when the 
energy of the impurity approaches the energy of the bottom of the subband. If the impurity 
energy decreases further, the current decreases but does not vanish. We argue that this 
phenomenon is due to phonon- assisted tunneling processes. Other additional features like 
Fermi edge singularities and fluctuations of the density of states are clearly resolved. We 
present an original analysis which unambiguously demonstrates that the peak at the Fermi 
level is related to the Fermi edge singularity and not to the fluctuations of the density of 
states. 
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FIG. 1: Sketch of the conduction band under applied bias. Solid line: F valley, dotted line: X 
valley. 



13 




0.5 1 1.5 2 2.5 3 



bias (V) 

FIG. 2: I{V) curves for large (0.5mm wide) doped (solid) and undoped (dashed) 
mesas. 
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FIG. 3: I{V) curve for a 400nm mesa measured at very low temperature 
(30mK). 
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FIG. 4: Detailed measurement of step 3 from Fig. 3. 
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FIG. 5: I{V) curves at T=30mK in the presence of a magnetic field of 0, 1, 1.5, 2 and 2.4 T. The 
field is oriented along the direction of the current. Curves have been vertically shifted, the shift 
is proportional to the magnetic field. Dotted lines are guides for eyes for the second Landau level 
and the bottom of the conduction band. 
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FIG. 6: Phonon-assisted tunneling rates for both deformation potential and 
piezoelectric interactions. 
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FIG. 7: Solid line: computed I{V) curve taking into account the elastic tun- 
neling contribution and both phonon-assisted contributions. Dashed line: a 

phenomenological broadening of the emitter density of states has been added 
in the calculation. 
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FIG. 8: First KV) step for both increasing and decreasing biases. 
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FIG. 9: Third I(y) step for both increasing and decreasing biases. 
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FIG. 10: /(T^) curves for both increasing (black line) and decreasing (gray line) biases, best 
obtained with FcSsCOlmeV. 
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